Negative Elliptic Flow of J/^'s: A Qualitative Signature for Charm Collectivity at 

RHIC 
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We discuss one of the most prominent features of the very recent preliminary elliptic flow data of 
J/ip meson from the PHENIX collaboration [J. Even within the the rather large error bars of the 
measured data a negative elliptic flow parameter (v?) for J/i/j in the range of pr = 0.5 — 2.5 GeV/c 
is visible. We argue that this negative elliptic flow at intermediate pr is a clear and qualitative 
signature for the collectivity of charm quarks produced in nucleus-nucleus reactions at RHIC. Within 
a parton recombination approach we show that a negative elliptic flow puts a lower limit on the 
collective transverse velocity of heavy quarks. The numerical value of the transverse flow velocity /3t 
for charm quarks that is necessary to reproduce the data is /3T(charm) ~ 0.55 — 0.6c and therefore 
compatible with the flow of light quarks. 



The main goal of the current and past heavy ion pro- 
grams is the search for a new state of matter called 
the Quark-Gluon-Plasma (QGP) Q- Major break- 
throughs for the potential discovery 0, 0] of this new 
state of matter were the observation of constituent quark 
number scaling of the elliptic flow v^ 100 - (pj? dron ) = 
UqV^iPx • on /n q ), with n q being the number of con- 
stituent quarks in the respective hadron as well as the 
observation of jet quenching at intermediate transverse 
momenta 0, @, 0|. Together with the 'standard' hydro- 
dynamical interpretation this implies a rapid thermal- 
ization and a strong collective flow of the QCD matter 
created at RHIC. However, open questions remain: how 
can one obtain a consistent description of the high px 
suppression and the elliptic flow of heavy flavour quarks 
and hadrons. I.e. is the collectivity at RHIC restricted to 
light quarks (up, down, strange) or do even charm (bot- 
tom) quarks participate in the collective expansion of the 
partonic system and reach local kinetic equilibrium? 

Previously, it was assumed that local equilibrium of 
(heavy) quarks could not be achieved within^ pQCD 
transport simulations. In fact, older studies [8| based 
on a parton cascade dynamics restricted to 2 <-> 2 parton 
interactions seemed to indicate that the opacity needed 
to achieve local equilibrium would be at least an order 
of magnitude higher than pQCD estimates. However, re- 
cent state-of-the-art parton cascade calculations (includ- 
ing 2 <-> 3 parton interactions) have clearly shown that 
pQCD cross sections are sufficient to reach local (gluon) 
equilibrium and allow to describe the measured elliptic 
flow data @, M, EI (3 • 

The aim of the present letter is to investigate whether 
also the charm quark does locally equilibrate and there- 
fore follows the flow of the light quarks. Here we will 
focus on the J/ij) because it reflects the momentum dis- 
tribution of the charm quarks directly, in addition first 
experimental data on the J/ip elliptic flow just became 
available. We will show that the recently measured nega- 
tive elliptic flow of J/ip's provides a unique lower bound 
on the charm quark's collective velocity. 

Under the assumption of local equilibration of light 



quarks a hydrodynamic parametrization of the freeze-out 
hyper-surface to parametrize the quark emission func- 
tion, namely the blast- wave model, can be employed. For 
the charm quarks, the same emission function is used, 
however, with the transverse collective velocity as a free 
parameter to be determined by the preliminary PHENIX 
data. To calculate J/^'s from the charm quark emission 
function, we a pply the well known parton recombination 
approach (13. Il4l fl5T |. Details (like the exact form of 
the freeze-out hyper-surface) of the specific approach em- 
ployed here can be found in [Ti.[l5|. Different from there 
we used a linear increasing transverse flow rapidity in- 
stead of a constant one, but the mean value has been 
preserved. 

Here we summarize the most important features: In a 
coalescence process the quarks contribute equally to the 
hadrons momentum, so it inherits its azimuthal asym- 
metry directly from its constituents. Therefore in recom- 
bination the elliptic flow of J/^'s emerges directly from 
a negative i>2 of the charm quark. To incorporate the 
asymmetry, the transverse expansion rapidity tjt depend 
on the azimuthal angle <j) and the radial coordinate p = 
as 

77t(0, p) = r,° T - \p (1 + ef(p T ) cos(20)) (1) 

with the eccentricity s and f(pr) = 1/(1+ (pt/po) 2 ) to 
model the damping at high px- With the factor | we re- 
cover rjj, as the mean transverse rapidity after integrating 
over p. 

By applying the definition of the elliptic flow one ob- 
tains fl3 | 
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f cos(20)/ 2 [a(<j),p)} Kx [6QM] d^palp 
J lo [a{<)>, p)} K x [b((/), p)] d(j> p dp 
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with a{(j>,p) = prsmh(r] T (^,p))/T, b{(j>, p) = 
rriT cosh(r)T(4>, p))/T and the modified Bessel functions 
I n and K n . For a more general hydrodynamical hyper- 
surface one could assume a dependence of freeze-out time 
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r on the radial coordinate p. This would lead to addi- 
tional terms involving |^ and Bessel functions of other 
order. We have checked that the modifications are only 
minor and therefore neglect the contributions in this let- 
ter for brevity. 

Let us investigate the elliptic flow of the J/tp at 
midrapidity as a function of the transverse momen- 
tum for various transverse flow velocities as shown in 
Fig. [TJ The lines from top to bottom indicate calcu- 
lations with a charm quark mass m c — 1.5 GeV/c 2 for 
different mean expansion velocities t&nh(r) T ) = f3 T = 
0.4c, 0.5c, 0.55c, 0.6c, the data by the PHENIX collab- 
oration are shown as symbols with error bars indicating 
a negative elliptic flow for J/ip's at intermediate trans- 
verse momenta. The calculation shows that with in- 
creasing transverse flow a negative vi at low pr (above 
Pt ~ 2.5 GeV,the elliptic flow values turn positive again) 
develops for the J/tp, posing a lower bound of fir > 0.5c 
for the charm quarks flow. The best fit to the data is 
obtained with a mean charm flow velocity of (5t = 0.55c 
equal to the light quark flow velocity extracted from pre- 
vious fits within the same model. 

In Fig. [2] we use /3r = 0.55c and compare the elliptic 
flow to other heavy mesons and Fig. [3] shows the same 
for the quarks. The value for T, with a bottom quark 
mass of m& = 4.7 GeV/c 2 , is negativ in the whole range 
of applicability. In contrast to J/tp , the t>2 of D° stays 
positiv. This is due to the positive light quark v%, which 
competes with the negative one for the charm quark, and 
results in nearly zero elliptic flow at low pr- While the 
B + meson follows the D° flow for pr < 1 GeV, it is much 
more suppressed at higher pr due to the strong negative 
flow of the bottom quark and approximately zero up to 
p T ~ 2.5 GeV. 

Data on £>-meson elliptic flow is not yet available. 
When comparing it to the non-photonic electron t>2, our 
calculations fail to predict the data flij |. These two ob- 
servables have been predicted to be similiar [13], since 
the non-photonic electrons are mainly from , -D-meson 
decays, but with a small contribution of S-meson de- 
cays. But the electron elliptic flow is no straightforward 
probe for the D° u 2 . Since the electron is not the only 
decay product, the decay kinematics might smear out 
the resulting elliptic flow of the electrons. At low pt, 
the increase of the D° flow is similiar to the electron 
data, but shifted to higher transverse momenta. Above 
Pt = 2 GeV the electron x>i starts to decrease which 
might be due to contributions from the i?-mesons or an 
early onset of the fragmentation regime. 

Direct measurements on the elliptic flow of heavy 
mesons will be available in the near future with the 
heavy-flavor tracker for STAR, which will allow a bet- 
ter analysis. Therefore the presented results are based 
only the J/tp elliptic flow data. 

These results provide strong evidence for a substantial 
collectivity and transverse expansion of the charm quarks 
in nucleus- nucleus reactions at RHIC. Due to the large 
error bars this has to be verified when more precise data 



is available. Note that our present findings are different 
from previous approaches that assume incomplete ther- 
malization of the charm [13, OH H3, H3] . We also verified 
our findings within a boltzmann approach to coalescence 
[20I ] using our parametrizations and received similar re- 
sults. 

One should also note that the observation of negative 
elliptic flow of heavy particles is well known in the litera- 
ture (even if not conclusively observed experimentally up 
to now). It appears due to an interplay between trans- 
verse expansion and particle mass, the more flow and the 
heavier the particle the more negative values does the el- 
liptic flow reach. E.g., negative values of the elliptic flow 
parameter for heavy hadrons has also been found in pre- 
vious exploratory studies and seem to be a general feature 
of the blast-wave like flow profile at high transverse ve- 
locities 0,0000]. lt reflects the depletion of the 
low pt particle abundance, when the source elements are 
highly boosted in the transverse direction. The difference 
to the present study is that here, vi is already negative 
on the quark level. Negative elliptic flow values will even 
be encountered for light quarks at asymptotically high 
bombarding energies as discussed in [lol |. One might ar- 
gue that this is an artefact of the blast-wave peak and 
will not survive in more realistic calculations, however 
also transport model calculations show slightly negative 
V2 values for heavy particles at low transverse momenta 

m 

In conclusion, we have shown that the recent prelimi- 
nary PHENIX data exhibiting a negative elliptic flow at 
low pt can be explained within a parton recombination 
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FIG. 1: Elliptic flow (v 2 ) of J/ip's for b = 9 fm for different 
mean transverse expansion velocities (lines) compared to pre- 
liminary data from PHENIX collaboration [J. While the V2 
of J ftp's is smaller than for light hadrons, the mean transverse 
velocity for the best-fit case ((3 = 0.55c for charm quarks) is 
the same as for light quarks. 
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FIG. 2: Comparison of elliptic flow (1)2) for J/^j £>°, T and 
B + at 6 = 9 fm with f) = 0.55c to data of non-photonic 
electrons from PHENIX collaboration (T3 |. 



approach using a blast-wave like parametrization. We 
point out that studying v 2 {pt) from J/tp offers the pos- 
sibility to put a lower limit on the charm quark trans- 
verse velocity. From the present quantitative analysis 
we expect the transverse velocity of charm quarks to be 
above (3t > 0.4c. Within the limits of the present model 
the best description of the data is obtain for a charm 
transverse velocity equal to the light quark velocity of 
/3r = 0.55c. So if more precise data will still support 
the negative v 2 , we conclude from this observation that 
charm quarks reach a substantial amount of local kinetic 
equilibration. 



0.15 - 



> 



0.05 - 



-0.05 - 



— I 1 1 

light quarks 

charm quark 
bottom quark 




_L 



1 2 3 4 5 

P T [GeV] 

FIG. 3: Elliptic flow (v%) of light, charm and bottom quarks 
at b = 9 fm with (3 = 0.55c. 
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